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Abstract 

Massive gravity (niGR) describes a dynamical “metric” on a fiducial, background one. We 
investigate fiuctuations of the dynamics about mGR solutions, that is about its “mean field 
theory”. Analyzing mean field massive gravity (mGR) propagation characteristics is not only 
equivalent to studying those of the full non-linear theory, but also in direct correspondence 
with earlier analyses of charged higher spin systems, the oldest example being the charged, 
massive spin 3/2 Rarita-Schwinger (RS) theory. The fiducial and mGR mean field back¬ 
ground metrics in the mGR model correspond to the RS Minkowski metric and external 
EM field. The common implications in both systems are that hyperbolicity holds only in 
a weak background-mean-field limit, immediately ruling both theories out as fundamental 
theories. Although both can still be considered as predictive effective models in the weak 
regime, their lower helicities exhibit superluminal behavior: lower helicity gravitons are su¬ 
perluminal as compared to photons propagating on either the fiducial or background metric. 
This “crystal-like” phenomenon of differing helicities having differing propagation speeds in 
both mGR and mGR is a peculiar feature of these models. 


1 Introduction 


Consistency is a powerful tool for studying field theories. Already classically, there are 
stringent conditions that are extremely difhcult to fulhll for systems with spin s > 1, the 
most important exception being (s = 2, m = 0) general relativity (GR). Key consistency 
requirements are 

(i) Correct degree of freedom (DoF) counts. 

(ii) Non-ghost kinetic terms. 

(iii) Predictability. 

(iv) (Sub)luminal propagation. 

Requirements (i) and (ii) are closely related (as are (iii) and (iv)). Models whose constraints 
do not single out the correct propagating DoF suffer from relatively ghost kinetic terms: 
the relevant example here is the sixth ghost excitation that plagues generic massive gravity 
(mCR) theories [1]. The discovery that a class of niGR models satished requirements (i) 
and (ii) generated a revival of interest in massive spin 2 theories [a El la El El [7] even though 
failure of the propagation requirements (iii) and (iv) were long known to bedevil higher spin 
theories iii. 

The predictability requirement is that initial data can be propagated to the future of 
spacetime hypersurfaces. In PDF terms, this means that the underlying equations must be 
hyperbolic [lO]. The hnal requirement, that signals cannot propagate faster than light, can 
be imposed once the hyperbolicity requirement is satished. The classic example of a model 
that obeys requirements (i) and (ii) as well as (iii) but only in a weak held region, is the 
charged, massive, s = 3/2 RS theory. Curiously enough, the propagation problems of this 
model were hrst discovered in a quantum setting by Johnson and Sudarshan [TT] who studied 
the model’s canonical held commutators (this is easy to understand in retrospect, because 
held commutators and propagators are directly related II2I)- The hrst detailed analysis of 
the model’s propagation characteristics was carried out by Velo and Zwanziger; our aim is 
to reproduce their RS results in mCR, so we quote their 1971 abstract verbatim [8]: 

The Rarita-Schwinger equation in an external electromagnetic potential is shown 
to be equivalent to a hyperbolic system of partial diherential equations supple¬ 
mented by initial conditions. The wave fronts of the classical solutions are calcu¬ 
lated and are found to propagate faster than light. Nevertheless, for sufficiently 
weak external potentials, a consistent quantum mechanics and quantum held 
theory may be established. These, however, violate the postulates of special 
relativity. 
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In previous works we and other authors have shown that similar conclusions hold for the 
full non-linear mGR models 0 El ESI EH 0 ES]- These investigations rely on the method 
of characteristics, which amounts to studying leading kinetic terms and is thus essentially 
equivalent to an analysis of linear fluctuations around a mean field background. Since this 
mean field massive gravity (mGR) fluctuation model depends both on a background and a 
fiducial metric, it is in direct correspondence with the charged RS model. Hence, without any 
computation at all, one can readily predict that: (a) mGR loses hyperbolicity in some strong 
held regime and (b) in the weak held hyperbolic regime where predictability is restored, 
lower helicity modes have propagation characteristics dihering from maximal helicity ±2; 
thus superluminality with respect to (luminal) photons is inevitable. Apart from conhrming 
earlier conclusions in a very simple setting, our results give a precise description of mGR’s 
ehective, weak held, regime. 


2 Massive Gravity 


At its genesis, the hrst known non-linear mGR model of [16] was originally formulated 
in terms of dynamical and hducial vierbeine and /”*. It took some forty years for 
researchers—independently in an ehective held theory-inspired metric formulation—to dis¬ 
cover that this model was one of a three-parameter family [2] that avoided the sixth, ghost-like 
excitation of [T]. The action describing these hducial mGR models is given bjfl 


5'mGR[e,C<;; /] = “ / ^mnrse^ { i ^ 




-F —e”e7" -1- f f + MTT 


The parameter Go governs a standard cosmological term; this is required to obtain the Fierz- 
Pauli (FP) linearized limit when both the hducial and mGR backgrounds are Minkowski. 
When both the hducial and mGR backgrounds are Einstein with cosmological constant A, 
the model’s parameters must obey ^ = m? (Go + Gi + /S2 + Gs) and the linearized theory 
is FP with mass m|p := m^(Gi + 2 G 2 + SGs)- 

Varying the model’s dynamical helds gives equations of motion 


( 1 ) 

where tm '■= Cmnrs [Goe^e^’e®-hGi^^'e''/^-1-G2e”/^/^+ G 3 /”/'’/*] • Also, the Einstein three-form 
is dehned by := \emnrs^"'R'^^ and i?™" := da;™"-f-is the Riemann curvature; V is 

'^Here d is the exterior derivative and the dynamical vierbiene and spin connection (e,uj), are one-forms. 
We suppress wedge products unless necessary for clarity. 
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the connection of a;™”. The forty equations above are subject to thirty constraints that are 
spelled out in detail in [7]. In particular, these include the covariant algebraic relation^ 

pm r n mrn^ 

^ Jm ~ ~ j.\.Ynn^ J i 

where the tensor := (jmn_^mn (jg^otes the contorsion and M™"" := /3ie™'e"'+2/32e[”*/"'] + 


3 Mean field massive gravity 


Consider mGR propagating in an arbitrary fiducial (pseudo-)Rieniannian manifold 
with corresponding vierbeine and spin connections (/”^, y™'"'). Now let (e™', cu™'"') be a solution 
to the mGR equations of motion ([1]). We wish to study fluctuations (e”^, A™'”) about this 
configuration: 

~m I ~,mn , ,mn i \ mn 

6 — 6“rc, UJ — UJ ~r A 

The action governing these is the quadratic part of S'mGR[e)i^; /] “ *S'mGR[e, ti;; /], namely 


S[h, A; e, /] := J e„ [e^^e^VA^^ + ^ 


The mean held model is a theory of forty dynamical helds (e™', A™"). In the above, V is the 
Levi-Civita connection of e”^, and its Riemann tensor; we stress that henceforth the 
hducial held (/”^, x”^"(/), ^/xi/(/)) and mGR background helds (e™, a;"*”(e), are non- 

dynamical] all index manipulations will be carried out using the mGR background metric 
and vierbein. 

The mGR equations of motion are 

;= + A™"e„ ^ 0 , 

Gm ■■= + ( 2 ) 


where := emnrs 


3/3o + 2/3i + ^2 


^The first of these assumes invertibility of the operator M"*" as a map from two-forms to antisymmetric 
Lorentz tensors; we shall always work in the model’s branch where this holds. 
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4 Mean field degrees of freedom 

In principle, since we are describing the linearization of a model whose constraints have been 
completely analyzed in [7] , we know a priori that mGR describes hve propagating degrees of 
freedom. However for completeness and our causality study, we reanalyze its constraints. 

The hrst step is to introduce a putative choice of time coordinate t, which for now need 
not rely in any way on either the hducial or background metric, and use this to decompose 
any p-form 0 (with p < 4) as 

e:=e + §, 

where 6 A dt = 0. Thus 6 is the purely spatial part of the form 0. Hence for any on-shell 
relation V ^ 0 polynomial in (V,£, A), its spatial 7^ ^ 0 part is a constraint because it 
contains no t-derivatives. 

Thus we immediately hnd sixteen primary constraints: 

r- = Ve’" + = ^emnrs - m . 

There are ten secondary constraints in total: The hrst six of these follow from the integrability 
condition ~ 0 which yields the so-called symmetry constraint 


T £[min] ~ 0 • 

As mentioned above, we assume that the set of two forms {M™”} is a basis for the space of 
two-forms, so the symmetry constraint yields 

e^fm ~ 0 . 


The remaining four secondary constraints come from the covariant curl VQm ~ 0 and give 
the vector constraint 

\7T -\- \ 0 

Employing the equation of motion ~ 0, this implies 

emnrs + 2(/3ie™ -F P2r) ~ 0 . 

Finally there are /our tertiary constraints stemming from covariant curls of the secondaries: 
the temporal part of the curled symmetry constraint + Xmn^'^f^ ~ 0, i.e. 


Jin I Ty- °m , ry rn i \ jin i \ °m jin , \ Jtr 

■^mn^ J -^mn^ J J \ J “r J J 
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and the scalar constraint 


emnrs - 2^28^^ f^) 

+emnrs - 2^26^^ - S/^s/™/*) 

+ 2 m^ + 2 + 2/32/"*) + 3 CmnM^rVX^^ 

-4/32£"*G™ - 26„A£'”e"i?*'^ ^ 0 . 

The VA'’® term seems to indicate that the above display is not a constraint for fits ^ 0, 
however as shown in [7], this quantity (weakly) equals one without time derivatives of helds. 
In summary, the model describes forty helds subject to thirty constraints and thus propagates 
hv^ DoF. 


5 Characteristic matrix 


We now study whether mGR can propagate initial data off a given hypersurface S. This 
amounts to asking if derivatives normal to S are determined by the equations of motion (|2]). 
For that, we simply replace all derivatives in the equations of motion and gradients of their 
constraints by the normal covector to S multiplied by the normal derivative of the corre¬ 
sponding held: 




= and 


a^A"*" 


= iadr^y 


We will also, for reasons of simplicity alone, restrict to the parameter choices /32 = /^s = 0 
(the model’s characteristic matrix for its entire parameter range has been computed in [7]). 
In particular we must focus on the question whether the linear system of equations for the 
normal derivatives := {d^e, c/nA) implied by the equations of motion along S is invertible. 
This amounts to a matrix problem encoded by the theory’s characteristic matrix C. In 
what follows we compute the system of equations given by the homogeneous linear system 
C ■ = 0. Starting with the equations of motion we hnd 


eAan£"* = 0, (3) 

emnrse"- A ^ A ^nA*"^ = 0 . 

^Although this conclusion for mGR is guaranteed by previous studies [HE] of the non-linear model’s 
DoF count, it verifies the linearized mGR study m- In that work, the fiducial metric is eliminated in terms 
of the mean field in order to argue that linearized spin 2 fields can propagate consistently in any gravitational 
background. This result is consistent with earlier work in |18j which relies on a 1/m^ expansion to study 
leading DoF and causality properties of gravitating, massive spin 2 models. 
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The gradients of the secondary and tertiary constraints then imply 


e/. fm A = 0 = emnrse^ ^ [e" A A 9^6*] , (4) 

/™ A [Kmn A + c” A dnXmn] = 0 , 

iymnrse^X {e ^+ {K^^AKi - + m2(4/3oe"Ae^ + 3(3ie^ Af^)) Adr,e^] = 0 . 

In the above the prefactor was included to indicate the origin of these equations but 
can be removed with impunity. To handle Equation ([3]) we decompose form-valued normal 
derivatives as earlier in Equation (|1]), and fincu 


Te™ = 0 = dr^y 


Supposing that the one-form ^ = dt, for some evolution coordinate t, we now use a shorthand 
notation = dte^ and = dtX^'^. We thus have the reduced characteristic system 


/ fm 

/ 9/^ V 

I ^^mnrs^ ^ 

1 /” X Knm 

\ 'Tim 




0 . 


(5) 


In the above square matrix, x denotes the standard three-dimensional cross product while 
the spatial densities on its last line can be read off from (jl]) and are simple for flat hducial 
metrics. Vanishing of the determinant of the above 10 x 10 matrix completely characterizes 
the boundary of the model’s predictive hyperbolic regime (modulo the restriction explained 
in Footnote HI). As we shall see, the reduced characteristic system describes the propagation 
of superluminal lower helicity modes: In the next section, we specialize to flat fiducial spaces 
and show how to analyze this determinant in direct analogy with the RS system. 


6 Analogy with Rarita—Schwinger 

The charged, spin 3/2, Rarita-Schwinger (RS) equation of motion reads 

771 

+ ieA^ + _ = 0 . 

Here V is the Levi-Civita connection of the fiducial spacetime and A is the background 
EM potential. These are analogous to the mGR hducial and background-mean-held metrics. 

^ Here we assumed that the pullback of the mGR background vierbeine to the hypersurface S is invertible 
and ignore likely pathologies when this fails. 
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The RS characteristic matrix was computed in [S] for flat hducial metrics and spacelike 
hypersurfaces S and found to have zero determinant when the magnetic held B obeyec@ 


-( 


2e y 




= 0 


( 6 ) 


The condition < 1 thus determines the weak held, hyperbolic, regime. Our aim now is 
to develop the analogous statement for mGR and inherit the conclusions of [8]. 

We begin with a short calculation. Consider now a hat hducial metric so 
and ds"^ = —dt^ + dx^. This simplihes the reduced characteristic system considerably. Firstly 
the equation fm^T = 0 of (E]) implies = 0, where we have decomposed the Lorentz index 
m = (0,a). Let us introduce an EM-like notation 


eabcK^^=:Ba, =: E\ 


For pure simplicity reasons only, we now restrict to the case where e° = E°' = 0 . Thus using 
e°abc = ^abc the secoud equation of ([5]) implies 

\ Oa ^ abc ri = ^0 

\ ~ ~ 2 ^ -Dft-ecG- 

Here the 3-vectors form the 3-inverse of e“ so that = 5^. The third equation of ([5]) 

then gives 

f[a X eb]Af = 0 . 

This equation generically allows to be expressed as a function of but this requires 
a non-trivial condition on e“. Under the hypothesis that the eigenvalue spectra of the 
matrices /“ ■ and —/“ ■ do not intersect, the above equation implies that A“^ = 0. In 
this framework, the last equation of ([5]) gives the single VZ-type condition 

mfp(4 - (r ■ e.)) - i((B“ ■ ■ e„) - i(B“ ■ ei.XB'’ • e,,)) 

In the weak held limit where the mean held approaches Minkowski space, the coefficient 
of rripp approaches unity but can change sign in a strong-held, large e limit. Hence there 
are certainly strong held conhgurations where the model loses hyperbolicity ra and closed 
causal curves are unavoidable. (This signals the onset of strong coupling in an ehective 
held theory.) Now, comparing Equations (E]) and ([7]), we see that we have reduced mGR’s 
weak held propagation analysis to a previous—well understood—case. Finally, we note that 
the same analogy and characteristic method can also be applied to the bimetric theory by 
treating the two background metrics as a (hducial,background) pair [20] . 

®See [19] for models designed to cure this pathology by adding higher background derivative, string- 
inspired terms to the RS action. 


e? = 0 . (7) 







7 Conclusions 


mGR is not a fundamental theory but rather an effective one with a range of validity de¬ 
termined by requiring hyperbolicity in a weak held regime. Excepting the further caveats 
explained in the text, the reduced characteristic matrix of Equation ([5]) completely deter¬ 
mines this allowed regime. Even in the weak regime, modes exhibit a crystal structure with 
differing maximal propagation speeds. For mGR to give a useful effective theory for physical 
applications, one must couple to matter (or at least photons) and require consistent causal 
cones for all modes. Once these couplings are decided upon, the characteristic method will 
determine their effective range of validity, if any. There is also the logical possibility that 
there exists a causal, luminal, UV completion of mGR analogous to that for QED in curved 
space [21] (see however the more general discussion of UV completions [22] ) • 
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